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We discuss a general fourth-order theory of gravity on the brane. In general, the formulation of 
the junction conditions (except for Euler characteristics such as Gauss-Bonnet term) leads to the 
higher powers of the delta function and requires regularization. We suggest the way to avoid such a 
problem by imposing the metric and its first derivative to be regular at the brane, while the second 
derivative to have a kink, the third derivative of the metric to have a step function discontinuity, and 
no sooner as the fourth derivative of the metric to give the delta function contribution to the field 
equations. Alternatively, we discuss the reduction of the fourth-order gravity to the second-order 
theory by introducing an extra tensor field. We formulate the appropriate junction conditions on 
the brane. We prove the equivalence of both theories. In particular, we prove the equivalence of 
the junction conditions with different assumptions related to the continuity of the metric along the 
brane. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Brane universes have made great popularity during the 
last years [1, 2]. However, it is remarkable that so far only 
the standard Einstein gravity, Gauss-Bonnet gravity [3- 
9] and, in general, Euler density gravity [10] on the brane 
have been considered in the literature [11]. These can be 
expressed by the general action 1 

s= [ d D xv^y2 Kn i {n) +s brane +s m , (i.i) 

where is the Euler density of the n-th order, n n is 
an appropriate constant of the n-th order, M is a D- 
dimensional manifold, Sbrane is the brane action and S m 
is the matter action. The lowest order Euler densities 
are: the cosmological constant = 1, the Ricci scalar 
= R, and the Gauss-Bonnet density 1^ = Rqb = 
R a bcdR abcd - 4:R ah R ab + R 2 with appropriate constants 
k = -2A(2k 2 )- x = -2A/16ttG, m = (2k 2 )" 1 , k 2 = 
a(2n 2 )~ 1 , a = const, etc., a, 6, c = 0, 1, . . . , D — 3, D — 
2,D [5]. 

In fact, in a general class of brane models based on 
an arbitrary combination of the higher-order curvature 
terms f(R abcd R abcd ,R ab R ab ,R) the field equations are 
fourth-order. Because of that they are plagued by the 
higher power terms of the second derivative of the warp 
factor function o~(y) =\ y |. This leads to a production 
of the higher powers of the delta function S(y) which can 
make the field equations ambiguous. Among the general 
class, the models based on the Euler densities are unique 
in the sense that the higher powers of the second deriva- 
tive of the warp factor d 2 a(y)/dy 2 exactly cancel in the 
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1 We use convention (-++...+) for the metric following Ref. [12]. 



field equations [11]. One then is easily able to formu- 
late appropriate junction conditions given first by Israel 
[13, 14]. 

The main objective of our paper will be the study of a 
general fourth-order theory of gravity on the brane [15] 

S = x- 1 / d D xV=gf(X,Y,Z) + S brane + S m (X2) 

JM 

where X = R, Y = R ab R ab , R abcd R abcd are curva- 
ture invariants, and \ is a constant. It includes the 
Euler density theories with the first Euler density be- 
ing just f(X, Y, Z) — x K iA = x K iR an d the second 
Euler density being the Gauss-Bonnet term given by 
f(X, Y, Z) = X K2{Z -4Y + X 2 ) etc. 

Up to our knowledge, the only non-Eulerian density 
cases were studied in Refs. [16] and [17]. In Ref. [16] 
the fourth-order theory f(X, Y, Z) = f(X) = f(R) was 
first reduced to the second-order theory, and then trans- 
formed into the Einstein theory. The junction conditions 
were then obtained, and they were obviously free from 
the problem of the powers of 5— function contribution. 
On the other hand, in Ref. [17] the theories with the lin- 
ear combination of the form f(X, Y, Z) = aX 2 + bY + cZ 
(a,b,c = const.) were considered and the junction con- 
ditions were obtained by the application of the appropri- 
ate Gibbons-Hawking boundary terms, again after trans- 
forming this theory to an equivalent second-order theory. 

It is important to emphasize that the theories based 
on the functions of the Euler densities such as f(I^) 
are fourth-order. Among them the most popular are 
f(I^) = f(R) - the theories of the function of the first 
Euler density [18]. In fact, the theories which are based 
on the function of the second Euler density f(I^) have 
also gained some interest recently [19], but they have not 
been studied on the brane yet. 

Our paper is organized as follows. In Section II we dis- 
cuss the main obstacle to formulate junction conditions 
for the fourth-order brancworld in a standard way which 
has been performed in the case of Euler densities. In Sec- 
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tion III we make a proposal to formulate these junction 
conditions by imposing more regularity onto the metric 
tensor. Since it does not necessarily satisfy everybody's 
taste we present in Section IV an alternative approach. 
In this approach we transform our general fourth-order 
theory into a second-order theory by applying generalized 
Lagrange-multiplicr approach [18, 19, 26]. This method 
was successful in obtaining the junction conditions in 
f(R) theory [16] and in / = aX 2 + bY + cZ theory [17]. 
In the Section V we formulate the junction conditions for 
the equivalent second-order theory. Finally, in Section VI 
we give our conclusions. 



II. PROBLEM FORMULATING ISRAEL 
JUNCTION CONDITIONS IN A 
FOURTH-ORDER BRANE WORLD 



In order to discuss the problem let us begin with the 
standard D-dimensional brane theory [1] whose action is 
just (1.1) with only m ^ 0, i.e., 

S = ^ j d D Xy/-gR + Sbrane + S m (II. 1) 

with the field equations 

G a b = R a b -(l/2)6 a b R=K 2 T a b , 
where the energy-momentum tensor is given by 
T b = T b -Q(-w) + T b +Q(w) + S(w)S a b , 



(II.2) 



(II.3) 



with S a b being the energy-momentum tensor on the 
brane, and T a b ± are the energy-momentum tensors on 
the both sides of the brane (i.e., in the bulk). 

For simplicity, let us assume that we work in Gaussian 
normal coordinates - i.e., that the D-dimensional metric 
is of the form ([i, v = 0, 1, 2, . . . , D — 2; w = D) 

ds 2 = g ab dx a dx b = edw 2 + h^dx^dx" , (11 A) 

where e = n-n — +1 for a spacelike hypersurface, e = — 1 
for a timelike hypersurface, and h a b = g a b — cn a rib is a 
projection tensor onto a (D — l)-dimensional hypersur- 
face, n is the normal vector to the hypersurface. The 
extrinsic curvature in these coordinates is defined as 



1 dh^ 



(11.5) 



2 dw 

By the application of the Gauss-Codazzi equations [23] 
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Rx^p = (D - 1) R\v, P + e[K lu ,K Xp -K ltp K Xv lT}..t 



one has the D-dimensional field equations in the form 

G w w = -\^R+ 1 -e[K 2 -Tr(K 2 )]= K 2 TUll.9) 

(11.10) 
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The Israel junction conditions can be obtained by the 
integration of the field equations (II.9)-(II.ll) in the limit 
lim„,^o f™ [3] and read as 
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(11.12) 
(11.13) 
(11.14) 



where [K^ v \ = + — . In general, for any quantity 
ft, one defines [ft] = fl + — Q~ , where fl^ means that this 
quantity was calculated on the left-hand-side and on the 
right-hand-side of the brane, respectively. 

The most important point is that these junction con- 
ditions are obtained provided we assume the following 
continuity conditions for the metric at w = [14]: 



h 



h 



= h+ v , (11.15) 
+ h+,, w , K~ v + K+ v , (11.16) 



which means that the metric is continuous at the brane 
but it has a kink, its first derivative has a step function 
discontinuity, and its second derivative gives the delta 
function contribution. In other words: 

VH - h- v (w)e(-w) + h+ v (w)6(w) , (11.17) 
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(11.18) 



(11.19) 



The equation (11.12) follows from (11.11) and (II. 3) as 
a consequence of the fact that the terms dK^/dw and 
dK/dw contain the delta function 6(w). However, in the 
fourth-order theory given by the action (1.2), the applica- 
tion of the continuity conditions (II.15)-(II.16) does not 
work. In order to discuss this let us first write down the 
field equations for the action (1.2) [15]: 



Pab 
pab 



X rp 

~e^-*-abi 

-\itf- 



fxR 



ab 



(11.20) 

2f Y R< a R b \ + 2f z R ed < a R b \ de 



+ fxMg ab g cd - 9 ac g hd ) + a(f Y R ab ) + g ab (f Y R c % d 



■- 2(.f Y R^ a ) b \-A(.f z R d ^ c ), cd , 



(11.21) 
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where fx = df/dX etc. The reason for not being the 
same continuity conditions (II.15)-(II.16) valid here is 
that the Riemann tensor 



R = V-VR+t 



dw 



where K = , Tr(K 2 ) = WK^, and appropriately, 
the Ricci tensor, and the Ricci scalar, as squared, give 
the terms 



0V" dK^ dK^ dK^ /dKY 
d 2 w dw ' dw dw ' \ dw J 



(11.22) 



which are proportional to S 2 (w), For example, in the sim- 
ple fourth-order gravity theory f(X, Y, Z) = f(X) = R 2 
the field equations contain the term fx,ww, which is pro- 
portional to S 2 (w). 

As mentioned in the Introduction, the only "mysteri- 
ous" cases which do not involve singular expressions are 
the Euler densities. In these cases the terms proportional 
to 8 2 (w) exactly cancel. These theories have been studied 
by a number of people (e.g. [3-8]) and their field equa- 
tions are given by varying the action (1.1) with m ^ 0, 
«4 7^ and k 2 ^ 



where 

Pppua Rppvcr ~t~ ^h^^R^p -\- 2/lp[ I ,_R (T ]y i 

+ Rhp W h a]p , (11.28) 

■V - ^KK^KZ + K ap K^K llv -2K lip KP-K av 



K 2 K^) . 



(11.29) 



Two different ways to obtain the junction conditions have 
been applied so far. In Rcf. [3] the limit lim^^o J_ w 
of the field equations (II.23)-(II.25) with the energy- 
momentum tensor given by (II. 3) has been taken. On the 
other hand, in Ref. [5] the appropriate Gibbons-Hawking 
[20] has been added to the Gauss-Bonnet action which 
cancelled normal derivatives of the metric variation. The 
form of the boundary term for this theory had been given 
earlier in Refs. [21, 22]. In fact, the exact form of the 
boundary term for the action (1.2) of this paper is not yet 
known and according to the claim of Ref. [21] it does not 
exist. According to Ref. [24] there is even an ambiguity 
of the choice of the boundary terms within Gauss-Bonnet 
brane world. A more general discussion of the boundary 
terms is given in Refs. [28, 29]. 

The discussion of the Gibbons-Hawking boundary 
term to the standard first-order Euler density (i.e. Ricci 
scalar) was also given, for example, in Ref. [25]. 
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The "dangerous" terms appear in the "correction" H a b 
to the Einstein equations which is 
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The junction conditions are obtained by using the conti- 
nuity conditions (II.15)-(II.16) and read [5] 

2a(3[J„„] - [J]V - 2[P] wva [KY a ) 
+ - [K] V - -^S^ , (11.27) 



III. A PROPOSAL TO FORMULATE 
CONSISTENT JUNCTIONS CONDITIONS IN 
THE FOURTH-ORDER THEORY ON THE 
BRANE 

Bearing in mind the roots of the irregularity of the field 
equations (II.20)-(II.21), and using the notation of [13], 
we suggest to define the following singular hypersurface 
of order three for which the metric derivatives and the 
extrinsic curvature on both sides of it are given by 
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i.e., the metric and its first derivative are regular, the sec- 
ond derivative of the metric is continuous, but possesses 
a kink, the third derivative of the metric has a step func- 
tion discontinuity, and as long as the fourth derivative of 
the metric on the brane produces the delta function con- 
tribution. On the other hand, according to [13] a singular 
hypersurface of order two is defined as 
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and it would describe the boundary surfaces character- 
ized by jumps in the energy-momentum tensor (e.g. the 
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FIG. 1: A schematic picture illustrating the domains of in- 
tegration used in derivation of the junction conditions. Here 
V = Gl + G2, Al = dGl + AO and A2 = 8G2 - AO. 



boundary surface separating a star from the surround- 
ing vacuum) . The physical interpretation of the singular 
hypersurfaces of order three is not so obvious, since it 
should be characterized by a jump of the first derivative 
of the energy-momentum tensor. 

In order to carry on let us notice that the field equa- 
tions (II.20)-(II.21) can be rewritten as 

ab Xrjiab/ 



V-gCa b w abd , d + ^gc ab v ab = ±T ab c ab ^ (in.8) 

where we have introduced is an arbitrary tensor field C ab , 
and 
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fxA9 ab g cd - 9 (ac g b)d ) + UyR a Y (ni.9) 

9 ab (f Y R cd );c - 2(f Y R d ^ - A{f z R d ^ c ), c , 



V ab = - 1 fg ab + f x R ab + 2f Y R c{a R b \ 



cde 



(111.10) 



contain third derivatives of the metric. In fact, W abd :d is 
proportional to S(w), and the energy-momentum tensor 
T ab is given by (II. 3). Now, we integrate both sides of 
(III. 8) over the volume V which contains the following 
parts (cf. Fig. 1): Gl, G2 - are the left-hand-side and 
the right-hand-side bulk volumes which are separated by 
the brane, Al = dGl + AO, A2 = 8G2 - AO are the 
boundaries of these volumes, and ^40 is the brane which 
orientation is given by the direction of the normal vector 
n. We have 
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Using the Gauss theorem first 
bd 



g(C ab W abd ). d dn= / ^C ab W aba n d da , 
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we then integrate in the limit V — > AO (which in Gaus- 
sian normal coordinates (II. 4) corresponds to the limit 
lim^o JZJ to obtain 



AO 



f=7C ab {[W] aM n d - |s Qb } = , (111.13) 



where [W} abd = W abd+ - W abd - . Since the tensor C ab is 
arbitrary one can then conclude that 



[W] abd n d - ^S ab = , 



(III. 14) 



and these are exactly the junction conditions for the the- 
ory under study. 

For example, in f(X, Y, Z) = f(X) = f(R) theory in 
D = 5 dimensions with metric 



ds 2 = -dt 2 (111.15) 
+ a 2 (t, w)[dr 2 + r 2 (d0 2 + sin 2 Qd4> 2 )] + dw 2 

they give 
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where (...)' = d/dw, a — a(w = 0) and the brane 
energy-momentum tensor = (—po,Po,Po,Po)- 



IV. EQUIVALENT SECOND-ORDER THEORY 
APPROACH 

It is known that the fourth-order gravity theory 

S = X- 1 J d D x^-gf(R) (IV A) 

is equivalent to a second-order theory with the action [18] 

S = X - 1 I d D x^g [f(Q) (Q-R) + f(Q)} (IV.2) 
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of which equation of motion is just Q = R, provided 
f"(Q) 7^ 0. In this approach f'(Q) may be interpreted 
as an extra scalar field (j> = f'(Q). 

Similar approach can also be used to the fourth-order 
gravity of the function of Gauss-Bonnet term 

S = j d D x^-gf{R GB ) (IV.3) 

by considering the second-order theory action [19] 

S = X- 1 J d D x^~g If (A) (R GB -A) + f(A)] (IV .4) 

The variation of this action with respect to A gives the 
equation of motion A = Rqb, provided f"(A) ^ and 
tp = f'(A) is interpreted as an extra scalar field. 

A step ahead can be made by considering more general 
fourth-order theory [26, 27] 

~9f{9ab,Rab) (IV.5) 



= J d D x^-g 



which can be made equivalent to the second-order theory 
provided that one introduces the tensor field 
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which reduces the appropriate field equations to the 
second-order and eventually transforms it into the Ein- 
stein theory with some new metric tensor. 

We follow these considerations by noticing that the 
theory given by the action (1.2) is a special case of a 
more general theory with the action of the form 
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In fact, it is enough to consider the theory (IV.7) in order 
to get all the previous actions (IV. 1), (IV.3) and (IV.5) 
as well as the theory with the linear combinations of cur- 
vature terms f(X, Y, Z) = aX 2 + bY + cZ, as in Rcf. 
[17]. 

In general, the action (IV.7) again gives immediately 
the fourth-order field equations, but there exist a possi- 
bility to formulate it in an equivalent way by using the 
action of the following form 



Si = X- 1 f d D 
Jm 

+ f{g a b,<i>cdef)} 
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where 4> abcd is the tensor field which is independent of 
g ab . The variation of the action (IV.8) with respect to 
4>abcd gives the equation of motion which is simply 

(IV.9) 
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The condition (IV. 10) holds for the fourth-order theory 
with the action (1.2), too. This means that (1.2) and 
(IV.8) are the equivalent actions. Introducing the tensor 

jjghij _ 9f(g a b 7 4>abcd) ^ 
d4>ghij 

and then varying the full action S = Si+S m with respect 
to both H abcd and g ab , we get the field equations in the 
following form: 
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where 
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V(g ab ,H cde f)= (IV.15) 

- ll" :l '>0 :ll „,l fl .,,.ll- h h + f(gab^klmn(gab,H cde f)) 

In fact, the possibility to express the fields (f> a bcd as a 
function of g ab and H cde f is guaranteed by the condition 
(IV.10). 



V. JUNCTION CONDITIONS AND THEIR 
CORRESPONDENCE IN THE EQUIVALENT 
SECOND-ORDER THEORY 

The tensor 

Mab 1 ab r , 
df 
= g J + A VH ( Pecd{gab-,tl ) 

2 og ab 

is a function of the fields g ab and H abcd only, and it does 
not give any contribution to the junction conditions. This 
is because H abcd are continuous on the brane, their first 
normal derivatives have a jump, and their second normal 
derivatives are proportional to the 5(y) function on the 
brane. Performing the same type of calculation as in 
the Section III the relation (IV.13) gives the junction 
conditions in the form: 
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Assuming that (IV.9) is fulfilled, and that 
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This means that the junction conditions (V.l) in the 
second-order theory (IV. 8) are completely equivalent to 
the juction conditions (III. 14) in the fourth-order theory 
(1.2). 

The equation (IV. 12), which corresponds to (IV. 9) im- 
plies that Ricmann tensor is continuous on the branc, 
since the fields 4> ab cd ar e regular functions of continuous 
fields H abcd and g ab . This, together with the Gauss- 
Codazzi equations (II. 6) imply the following junction 
conditions: 



[K ab ] = 
[C n K ab \ = , 
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where £3 is a Lie derivative of H a ( e ^> d along the vector 
field n. (Vector field n is a normalized vector field tangent 
to the Gaussian normal coordinate w). The equations 
(V.4) reconstruct the assumptions (V.l). 

Bearing in mind the formulas above, equations (V.l) 
can be reduced to: 

n.n^CnA^} ^ -^S ab , (V.5) 

where 

A abcd = f x (g ad g cb - g cd g ba ) 

+ f Y (2R ad g bc - R cd g ba - R ba g cd ) 

+ Af z R achd . (V.6) 

In Gaussian normal coordinates the conditions (V.4) 
and (V.5) read as 
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In a specific example of f(X, Y, Z) = Z theory for the 
metric (III. 15), the equations (V.4) and (V.5) give: 



Pa = -Po 



(V.9) 
(V.10) 



VI. SUMMARY 

In this paper we have considered the junction condi- 
tions for the fourth-order brane world given by the action 
which is a general function of the curvature invariants 
R, R ab R ab , and R abcd R abcd . We imposed the regularity 



conditions on the metric tensor which was taken to be of 
the class C 2 functions of the coordinates. Explicitly, it 
means that the metric and its first derivative are regu- 
lar at the brane, while the second derivative has a kink, 
the third derivative of the metric has a step function dis- 
continuity, and the fourth derivative of the metric gives 
the delta function contribution to the field equations. In 
terms of the seminal notation given first by Israel [13], 
these conditions describe the singular hypersurfaces of 
order three. The junction conditions which we obtained 
are quite generic and they may be applied to some special 
cosmological framework of interest. 

As an alternative which allows less restrictive regular- 
ity conditions, we considered the reduction of the fourth- 
order theory to a second-order theory by applying an 
extra tensor field. We then formulated the junction con- 
ditions within such a theory and showed that they were 
equivalent to the previously obtained fourth-order theory 
junction conditions. 

In the previously considered cases in the literature 
mainly the Eulerian density brane worlds were studied 
[3-9, 11]. The only non-Eulerian density cases were in- 
vestigated in Refs. [16] and [17], where the theories with 
at most a linear combination the curvature invariants of 
the form f(X, Y, Z) = aX 2 + bY + cZ (a, b,c = const.) 
were considered. In these references the junction condi- 
tions were obtained after transforming such theories into 
the equivalent second-order theories. In fact, we made 
one step further, by suggesting junction conditions for 
the brane world which allows a general function of curva- 
ture invariants f(X, Y, Z). Also, in one of our approaches 
to the problem of junction conditions we do not reduce 
the fourth-order action into the second-order action, but 
suggest junction conditions to work in the fourth-order 
theory, though in the case of the so-called singular hy- 
persurfaces of order three only [13]. 

We hope that our calculations will allow us to study 
some cosmological applications of the general fourth- 
order gravity on the brane (1.2) in the following papers. 



VII. ACKNOWLEDGEMENTS 

We thank Andrei Barvinsky, Salvatore Capozzicllo, 
Tomasz Dcnkicwicz, Krzysztof Meissner and David 
Wands for discussions. M.P.D. acknowledges the support 
of the Polish Ministry of Science and Higher Education 
grant No 1 P03B 043 29 (years 2005-2007). 



[1] L. Randall and R. Sundrum, Phys. Rev. Lett., 83, 3370 
(1999); L. Randall and R. Sundrum, ibidem, 83, 4690 
(1999). 

[2] M. Visser, Phys. Lett. B159, 22 (1985); N. Arkani- 
Hamed, S. Dimopoulos, and G. Dvali, Phys. Lett. B516, 



70 (1998); I. Antoniadis, N. Arkani-Hamed, S. Dimopou- 
los, G. Dvali, Phys. Lett. B436, 257 (1998); N. Arkani- 
Hamed, S. Dimopoulos, and G. Dvali, Phys. Rev. D59, 
086004 (1999); P. Binetruy, C. Deffayet and D. Lan- 
glois, Nucl. Phys. B565, 269 (2000); P. Binetruy, C. Def- 



7 



fayet and D. Langlois, Phys. Lett. B477, 285 (2000); 
M. Sasaki, T. Shiromizu and K. Maeda, Phys. Rev. 
D62, 024008 (2000); T. Shiromizu, K. Maeda, and M. 
Sasaki, Phys. Rev. D62, 024012 (2000); S. Mukhoyama, 
T. Shiromizu and K. Maeda, Phys. Rev. D62, 024028 
(2000); A.N. Alicv and A.E. Gumrukgiiglu, Class. Quan- 
tum Grav. 21, 5081 (2004); D. Cremades, L.E. Ibanez 
and F. Marchesano, Nucl. Phys. B 643, 93 (2002); C. 
Kokorelis, Nucl. Phys. B 677, 115 (2004). 
[3] N. Deruelle and T. Dolezel, Phys. Rev. D62, 103502 
(2000). 

[4] C. Charmousis, J.F. Dufaux, Class. Quantum Grav. 19, 

4671 (2002). 
[5] S.C. Davis, Phys. Rev. D 67, 024030 (2003). 
[6] J.F. Dufaux, J.E. Lidsey, R. Maartens, and M. Sami, 

Phys. Rev. D 70, 083525 (2004). 
[7] J.E. Lidsey, Ann. Phys. (Leipzig) 15, 277 (2006). 
[8] H. Maeda, V. Sahni, Yu. Shtanov, gr-qc/0708.3287. 
[9] P.S. Apostopoulos et al, hep-th/0708.0469. 
[10] D. Lovelock, J. Math. Phys. 12, 498 (1971). 
[11] K.A. Meissner and M. Olechowski, Phys. Rev. Lett. 86, 

3708 (2001). 

[12] S.W. Hawking, G.F.R. Ellis, The Large-scale Structure 
of Space-time, (Cambridge Univ. Press, 1999) . 

[13] W. Israel, Nuovo Cimento B 44, 1 (1966). 

[14] M. Visser, Lorentzian Wormholes (Springer- Verlag, 
1996). 

[15] T. Clifton and J.D. Barrow, Phys. Rev. D 72, 123003 
(2005); T. Clifton and J.D. Barrow, Class. Quantum 
Grav. 23, 2951 (2006). 

[16] M. Parry, S. Pichler, and D. Deeg, JCAP 0504, 014 
(2005). 

[17] S. Nojiri and S.D. Odintsov, JHEP 0007, 049 (2000); S. 
Nojiri, S.D. Odintsov and S. Ogushi, Phys. Rev. D 65, 
023521 (2001). 

[18] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980); G. Mag- 
nano and L.M. Sokolowski, Phys. Rev. D 50, 5039 (1994); 



T.P. Sotiriou, Class. Quantum Grav. 23, 5117 (2006); 
T. Chiba, T.L. Smith, and A.L. Erickcek, Phys. Rev. 
D 75, 043516 (2007); G.J. Olmo, Phys. Rev. Lett. 98, 
061101 (2007); G.J. Olmo, Phys. Rev. D 75, 023511 
(2007); S. Capozziello, V.F. Cardone, and A. Troisi, 
Phys. Rev. D 71, 043503 (2005); S. Capozziello, S. NO- 
jiri, S.D. Odintsov, and A. Troisi, Phys. Lett. B 639, 135 
(2006); T. Chiba, Phys. Rev. D 75, 043516 (2007); S. 
Capozziello and R. Garatini, Class. Quantum Grav. 24, 
1627 (2007); L. Amendola, D. Polarski, and S. Tsujikawa, 
Phys. Rev. Lett. 98, 131302 (2007); T. Faulkner et al, 
astro-ph/0612569. 

[19] S. Nojiri, S.D. Odintsov, Phys. Lett. B 631, 1 (2005); S. 
Nojiri, S.D. Odintsov, O.G. Gorbunova, J. Phys. A 39, 
6627 (2006); G. Cognola et al, Phys. Rev. D 73, 084007 
(2006); B. Li, J.D. Barrow, D.F. Mota, Phys. Rev. D 76, 
044027 (2007); S.C. Davis, hep-th/0709.4453. 

[20] G.W. Gibbons and S.W. Hawking, Phys. Rev. D 15, 2752 
(1977). 

[21] T.S. Bunch, Journ. Phys. A 14, L139 (1981). 

[22] F. Miiller-Hoissen, Phys. Lett. B 163, 106 (1985); R.C. 
Myers, Phys. Rev. D 36, 392 (1987). 

[23] P.D. Mannheim, Brane Localized Gravity (World Scien- 
tific, 2005). 

[24] S. Nojiri and S.D. Odintsov, Gen. Rel. Grav. 37, 1419 

(2005) . 

[25] E. Papantopoulos, Lect. Notes Phys. 592, 458 (2002). 
[26] A. Jakubiec and J. Kijowski, Phys. Rev. D 37, 1406 
(1988). 

[27] J.D. Barrow and S. Hervik, Phys. Rev. D 73, 023007 

(2006) ; ibidem, 74, 124017 (2006). 

[28] A.D. Barvinsky and S.N. Solodukhin, Nucl. Phys. B 479, 
305 (1996). 

[29] T.P. Sotiriou and S. Liberati, Phys. Rev. D 74, 044016 
(2006). 



